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1.  INTRODUCTION 


The  primary  source  of  NO  emissions  in  an  aurorally  enhanced  upper  atmosphere  is  through 
odd-nitrogen  chemiluminescent  reactions  involving  N('*S),  N(^D),  and  N(^P).  Recent  papers 
speculating  on  the  creation  of  the  NO  hot  bands  have  identified  high  energy  nitrogen  atoms  as 
the  source  of  these  emissions.  One  of  the  goals  of  the  EXCEDE  III  measurement  program, 
which  has  been  described  previously, is  to  understand  the  electron-beam  induced  infrared 
fluorescence  resulting  from  chemical  reactions  and  energy  transfer  processes.  To  characterize 
the  relative  importance  of  the  N  atom  precursors  in  NO  formation,  the  N(^D)  and  N(^P) 
emissions  have  been  measured  in  EXCEDE  III.  These  measurements  are  compared  to  a 
chemical  kinetics  model  which  is  tested  for  consistency  by  using  supporting  measurements,  such 
as  electron  deposition,  NO  chemiluminescence,  and  N2  second  positive  (2P)  emissions.  The 
atomic  nitrogen  formed  when  an  energetic  electron  interacts  with  the  atmosphere  determines  not 
only  how  much  NO  is  formed,  but  the  spectral  distribution  of  the  chemiluminescence.  The 
relative  production  of  these  atoms  in  a  electron  excited  atmosphere  has  been  measured  under 
controlled  conditions  in  EXCEDE  III.  Simultaneous  measurements  of  the  ^P  ^  ^S,  ^D  ^  '^S,  and 
N2  first  negative  (IN)  bands  have  been  performed  during  this  experiment.  A  carefully 
controlled  electron  beam  dosing  profile  gave  a  direct  measure  of  the  odd  nitrogen  production 
rates  as  a  function  of  electron  injection.  The  energy  deposition  profile  was  measured  using  3914 
A  emission,  and  is  well  represented  by  a  generalized  Gaussian  function.  Since  the  exact 
mission  profile  was  known,  the  emissions  within  the  field-of-view  of  each  detector  could  be 
compared  with  the  chemical  kinetics  responsible  for  the  nitrogen  formation.  This  report  presents 
an  analysis  of  the  nitrogen  atoms  produced  in  its  two  lowest  metastable  excited  states  and 
predicts  the  amount  of  N('*S)  formed  based  on  branching  ratios  reported  in  the  literature. 

The  instruments  used  to  obtain  the  N  atom  concentrations  consisted  of  a  visible  spectrometer 
for  N(^P),  and  dedicated  photometer  for  N(^D).  A  second  photometer,  looking  at  Nj  IN,  was 
boresighted  with  each  of  these  instruments.  Since  several  different  instruments  were  used  for 
these  measurements,  they  are  described  separately  in  the  sections  below.  In  the  following 
sections,  we  present  the  measurement  of  N(^D)  and  N(^P),  the  modeled  concentrations  based  on 
the  observed  Nj  (0,0)  IN  emissions,  and  finally  a  discussion  on  the  comparison. 
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2.  MEASUREMENT  OF  METASTABLE  N(^D)  AND  NC^P) 


2.1  Nr^D) 

The  N(^D)  emission  was  measured  by  a  filter  photometer  with  a  center  wavelength  at  5200 
A  and  a  full  width  at  half  maximum  (FWHM)  of  6.06  A.  The  field-of-view  (FOV)  of  the 
photometer  was  nominally  6°  and  was  pointed  18°  into  the  afterglow.  Data  were  collected  at 
a  sample  rate  of  750  Hz.  In  order  to  correct  for  possible  interfering  emission,  such  as  the  0-3 
band  of  the  N2  IN  series  at  5228  A  or  the  Herzberg  bands  of  O2,  an  additional  photometer  with 
a  center  wavelength  at  5228  A  with  a  FWHM  of  5.79  A  was  co-aligned  with  the  5200  A 
photometer. 

2.2  N(¥) 

The  measurement  of  N(^P)  produced  in  EXCEDE  III  was  made  by  measuring  its  relative 
intensity  to  known  features  using  a  visible  spectrometer.  A  more  detailed  description  of  this 
instrument  is  reported  by  Reider  et  al.^^^  This  instrument  had  a  spectral  scan  time  of  2.8  s. 
Since  alternating  scans  used  a  neutral  density  filter  (OD=2)  to  extend  the  instrument  dynamic 
range,  one  3466  A  measurement  was  made  every  5.6  s.  This  time,  convolved  with  beam-on 
times  of  7.1  s,  limited  the  number  of  N(^P)  scans  to  15,  6  on  upleg,  one  at  apogee,  and  8  on 
the  downleg.  The  feature  at  3466  A  is  in  the  region  of  the  spectrum  where  the  Nj  and  N2  2nd 
Positive  bands  are  located.  In  fact,  the  Av  =  1  series  of  N2(2+)  has  a  band  in  near  coincidence 
with  N(2p),  at  3468  A. 
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3.  CHEMICAL  KINETICS  MODEL 


The  most  important  infrared  emitters  in  an  electron-dosed  atmosphere  include  NO,  CO2, 
NO'*',  CO,  and  O3.  A  general  chemical  kinetics  model  has  been  developed  to  describe  the 
chemical  and  collisional  processes  occurring  during  the  dosing  of  the  atmosphere  with  an 
energetic  electron  beam.  The  model  uses  a  finite  rate  approach  to  the  modeling  of  the 
generation  of  energy  dependent  secondary  electrons,  vibrationally  and  electronically  excited 
species  resulting  from  collisional  energy  transfer  and  chemiluminescent  reactions.  In  particular, 
time  and  energy  dependent  rate  equations  are  solved  for  arbitrary  electron  beam  dosing  to  obtain 
species  number  densities  as  a  function  of  time.  The  secondary  electron  spectrum  is  treated  by 
setting  up  an  energy  grid  (histogram)  at  sufficiently  high  resolution  to  accurately  solve 
electron-molecule  and  electron-electron  interactions.  Time-dependent  rate  equations  are  then 
solved  for  each  energy  bin.  Chemical  and  collisional  processes  for  individual  vibrational  states 
are  treated  explicitly.  Self-relaxation  of  secondary  electrons,  which  is  important  due  to  the  high 
electron  beam  fluxes  produced  in  the  EXCEDE  experiment,  has  been  treated  by  solving  the 
time-dependent  Fokker-Planck  equation. In  the  remainder  of  this  section,  the  detailed 
chemical  and  collisional  processes  influencing  the  production  and  destruction  of  metastable 
nitrogen  atoms  is  discussed  in  some  detail. 

When  electrons  interact  with  the  natural  atmosphere,  dissociation  and  ionization  form  odd 
nitrogen  in  the  form  of  N(^P),  N(^D),  and  Nf'^S).  Collisions  between  these  species  and  O2 
create  This  production  rate  is  in  turn  found  to  be  dependent  on  the  ratio  of  N(^D)  to 

N('*S)  formed  in  the  interaction  of  electrons  and  N2.  Nitric  oxide  is  currently  thought  to  be 
produced  primarily  by  the  reaction  of  N(^D)  with  O2 

NC^D)  -t-  O2  NO  -I-  O  (1) 

where  the  reaction  rate  constant^^^  and  NO  vibrational  distribution^^ are  well  established  at 
room  temperature.  Recent  laboratory  studies^®^  have  indicated  that  N(^P)  produces  NO  in  high 
vibrational  and  rotational  states 

N(2p)  -I-  O2  -►  NO  -I-  O  (2) 

with  a  rate  constant  a  factor  of  2  smaller  than  (1).^^^  Furthermore,  it  has  been  suggested  by 
Sharma  et  al.^^^  that  translationally  hot  NC'^S)  atoms^*®^  may  be  a  significant  source  of  rotationally 
hot  NO  in  the  thermosphere 
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N(^S)*  +  O2  -►  NO  +  O  .  (3) 

Although  the  reaction  rate  constant  and  the  final  vibrational-rotational  distribution  at  high 
translational  energies  have  not  been  measured,  recent  calculations  by  Duff  et  al/*^^  have 
established  realistic  rate  constants  and  NO  vibrational-rotational  distributions  for  Reaction  (3). 
A  realistic  description  of  the  formation  and  emission  of  NO  in  the  ambient  or  auroral  atmosphere 
requires  a  detailed  treatment  of  nitrogen  atom  chemistiy. 

Nitrogen  atoms  are  produced  by  several  major  processes;  dissociation  of  N2  by  the  primary 
electrons  (Cp)  in  the  electron  beam,  Reaction  (4),  dissociation  of  N2  by  the  secondary  electrons 
(Cj)  produced  by  ionization  of  atmospheric  species  by  the  electron  beam,  Reaction  (5),  and 
recombination  reactions  of  the  primary  ionic  products  with  secondary  electrons.  Reactions  (6) 
and  (7), 


ep  +  N2 

ep  -b  N(^S)  N(2D) 

(4a) 

Cp  +  N2 

ep  N(^S)  +  N(2p) 

(4b) 

Cp  +  N2 

ep  +  e3  +  N'*'  -b  ‘/2N(^S)  -b  ‘/2N(2D) 

(4c) 

Cs  +  N2 

e;  -b  N(^S)  -h  N(2D) 

(5a) 

es  +  N2 

e;  -b  N(^S)  -h  N(¥) 

(5b) 

NJ  +  63 

e3  -H  (2-ai)N('*S)  -H  aiN(2D) 

(6) 

NO+  +  63  -► 

(l-a2)N(^S)  -1-  a2N(2D)  -1-  0  . 

(7) 

The  total  production  rate  of  N  atoms  from  Reactions  (4a)  and  (4b)  is  assumed  to  be 
1.5P(N2),^^^’  where  P(N2)  is  the  Nj  production  rate,  with  the  relative  formation  rate  of 
N(^D)  and  N(^P)  in  the  ratio  0.52:0.48/^^^  For  dissociative  ionization  by  primary  electrons. 
Reaction  (4c),  the  production  rate  is  0.50P(N  Thus,  the  total  production  rate  of  N  atoms 

from  Reaction  (4)  is  2P(N2),  or  1.24  N  atoms/ion  pair,  in  good  agreement  with  the  results  of 
Porter  et  al.^^^^  and  Rusch  et  al.^^^^  The  cross  sections  for  dissociation  of  N2  by  low  energy 
secondary  electrons.  Reaction  (5),  are  based  on  the  analysis  of  Zipf  et  al.^’"^^  The  rate  constants 
for  the  recombination  of  Nj  and  NO'*'  ions  with  secondary  electrons.  Reactions  (6)  and  (7), 
have  been  measured  by  Mehr  and  Biondi^^^^  and  Alge  et  al.,^^^^  respectively.  The  associated 
branching  ratios  for  the  N2  and  NO'*'  recombination  reactions  are  aj  =  1.86(^*)  and  a2=0.76.(^^^ 
In  addition  to  Reactions  (4)  through  (7),  minor  sources  of  N(^D)  considered  are  the  charge 
exchange  reactions  of  Nj  with  and  N^  with  02^^*’^^^ 
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Nj  +  O  N(2D)  +  NO+  (8) 

N+  +  O2  N(2D)  +  OJ  .  (9) 

Nitrogen  atoms  are  also  subject  to  quenching  by  other  species  which  compete  with  Reactions 
(1)  and  (2)  in  the  production  of  NO 

N(2D)  +  N2  Nf'^S)  +  N2  (10) 

N(2D)  +  O  -►  N(^S)  +0  (11) 

N(2p)  +  O  N(2D)  +  O  (12a) 

N(2p)  +  0  NO+  +  Cs  .  (12b) 


The  quenching  rates  of  Reactions  (10)-(12)^^’^^^  are  slow  on  the  time  scale  of  the  EXCEDE  III 
measurements  in  comparison  to  the  efficient  quenching  of  N(^D)  and  N(^P)  by  secondary 
electrons  which  results  from  the  large  density  of  secondary  electrons  created  by  the 


ionization  rates  achieved  in  the  experiment 

N(2D)  +  Cj  -►  N(^S)  +  e;  (13) 

N(¥)  +  Ni%^D)  +  e;  .  (14) 

In  addition  to  collisional  quenching,  N(^D)  and  N(^P)  can  undergo  radiative  relaxation 

N(2d)  -►  N(^S)  +  hv  (15) 

N(2p)  -4  N(2d,^S)  +  hi'  .  (16) 


Since  the  N  atom  production  from  direct  electron  impact  excitation  and  recombination 
reactions  depends  on  the  low  energy  secondary  electron  spectra,  it  is  important  that  the 
histogram  treatment  retain  the  major  features  of  the  secondary  electron  spectra.  Secondary 
electrons  are  formed  via  collisional  ionization  of  M(=  N2,  O2,  or  O) 

Cp  +  M  Cp  +  Cg  +  products  .  (17) 

The  nascent  secondary  electron  energy  distribution  resulting  from  (17)  is  that  given  by  Opal  et 
al.,^^^^  which  then  relaxes  by  collisions  with  major  species 

e,  +  M  -►  e;  +  M*  (18) 

and 

e^  +  M  e;  +  e;'  +  M+  (19) 

as  well  as  electron-electron  collisions 

Cs  +  e;  e;'  -H  e”'  (20) 

where  the  primes  refer  to  secondary  electrons  in  different  energy  bins.  The  ion  production  rates 

from  Reaction  (17)  are  related  to  the  ion  pair  production  rate  via  the  expressions  given  by  Rees 
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and  Jones/^^^  The  rate  constants  for  Reactions  (18)  and  (19)  were  taken  from  the  reviews  of 
Itikawa  et  for  N2  and  O2  and  Itikawa  and  Ichimura^^^^  for  O. 

One  of  the  most  important  inputs  to  the  model  is  the  electron  beam  characteristics,  i.e.,  the 
radial  profile  of  the  ion-pair  production  rate.  Historically,  the  electron  deposition  information 
has  been  deduced  from  the  Nj  first  negative  3914  A  emission,  where  it  is  assumed  that  there 
are  14.1  Nj  ions  created  for  each  3914  A  photon.^^^^  Therefore,  the  beam  ion  pair  production 
rate  profiles  used  in  the  current  model  are  based  on  an  analysis  of  the  3914  A  data  by  Rieder 
et  al.^^^  The  Nj  production  rate,  P(N^),  is  related  to  the  total  ion  pair  production  rate,  Pio„, 
by 

P(NJ)  = 

1.15[N2]  +  1.5[02]+0.56[0] 

where  the  numerical  factors  in  the  expression  represent  the  relative  ionization  cross  sections  for 
N2,  O2,  and  Using  the  atmospheric  profile  appropriate  for  the  EXCEDE  III  experiment, 
the  relative  production  rate  of  Nj  ions  is  approximately  independent  of  altitude  (ranging  from 
0.60  at  95  km  to  0.62  at  115  km). 

Representative  predictions  of  the  production  and  destruction  rates  for  N(^D)  and  N(^P)  are 
presented  in  Figures  1  through  4  for  an  altitude  of  115  km.  The  electron  beam  ion  pair 
production  rate  is  represented  by  a  generalized  Gaussian  function^^^  with  a  peak  ion  pair 
production  rate  of  1.6xl0‘*  ion  pairs/cm^-s  and  a  full  width  at  half  maximum  of  14  ms. 
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TIME  (s) 

Figure  1.  Major  Processes  for  the  Production  Rate  of  N(^D)  at  115  km. 


TIME  (s) 

Figure  2.  Major  Processes  for  the  Destruction  Rate  of  N(^D)  at  115  km 


DESTRUCTION  RATE  (cm'^s'^l  ^  PRODUCTION  RATE  (cm'^s'^) 


4.  COMPARISON  OF  DATA  WITH  THE  KINETICS  MODEL 


The  chemical  kinetics  mechanism  described  previously  has  been  incorporated  into  the 
EXCEDE  chemical  kinetics  data  base.  The  resultant  data  base  contains  6703  reactions  and 
energy  transfer  processes.  The  chemical  kinetics  equations  are  integrated  using  Gear’s  method 
for  stiff  differential  equations^^®^  as  a  function  of  time  for  a  specified  electron  beam  ion-pair 
production  rate  and  altitude.  The  atmospheric  profiles  for  N2,  O2,  and  O  were  obtained  from 
the  MSISE-90  atmospheric  model^^^^  for  the  geophysical  parameters  appropriate  to  the  EXCEDE 
III  experiment. 

The  geometry  of  the  flight  experiment,  which  includes  the  trajectory  of  the  gun  module,  the 
geomagnetic  azimuth  and  dip  angles,  and  the  viewing  geometry  of  the  3466  A  visible 
spectrometer  and  5200  A  photometer  on  the  sensor  module,  has  also  been  incorporated  in  the 
model.  The  time-dependent  nitrogen  atom  number  densities  are  integrated  along  the  instrument 
line-of-sight  to  obtain  the  column  densities  as  a  function  of  time.  The  resultant  column  densities 
are  then  integrated  over  the  6“  FOV  of  the  instruments  for  comparisons  with  the  column 
densities  obtained  from  3466  A  and  5200  A  data. 

To  illustrate  the  electron  beam  energy  deposition  throughout  the  EXCEDE  flight,  the  total 
number  of  ion  pairs  produced  by  the  electron  beam  is  shown  as  a  function  of  mission  elapsed 
time  (or  altitude)  in  Figure  5.  The  number  of  ion  pairs  was  obtained  by  integrating  the  observed 
3914  A  scanning  photometer  data  along  the  component  of  gun  vehicle  velocity  vector 
perpendicular  to  the  magnetic  field  lines.  The  maximum  number  of  ion  pairs  are  created  near 
102  km  on  upleg,  where,  due  to  the  near  zero  velocity  across  the  field  lines,  the  atmosphere  is 
dosed  for  nearly  0.5  seconds.  The  average  electron-deposition  time,  which  is  simply  the  root 
mean  square  electron  beam  diameter  divided  by  the  velocity  across  the  field  lines,  is  shown  in 
Figure  6.  As  altitude  increases,  the  number  of  ion  pairs  produced  decreases  due  to  the  shorter 
deposition  times  (i.e.,  the  velocity  across  the  field  lines  is  increasing  with  altitude),  see  Figure 
6.  On  downleg,  the  number  of  ion  pairs  again  increases  as  altitude  decreases,  although  the  total 
number  of  ion  pairs  created  is  substantially  less  on  downleg  than  on  upleg.  This  effect  is  due 
to  the  atmosphere  being  irradiated  by  the  electron  beam  for  times  up  to  a  factor  of  50  longer  on 
upleg  than  on  downleg,  again  as  shown  in  Figure  6.  The  two  data  points  at  approximately  114 
and  113  km  are  about  75%  of  the  expected  beam  power  output  due  to  beam  load  faults. 
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The  comparisons  of  the  N(^D)  and  N(^P)  column  densities  obtained  from  the  5200  A  and 
3466  A  data  with  the  model  predictions  which  have  been  integrated  along  the  line-of-sight  and 
integrated  over  the  sensor  field-of-view  are  shown  in  Figures  7  and  8,  respectively. 

The  agreement  between  the  model  predictions  and  data  is  excellent  on  both  upleg  and  downleg 
for  the  N(^D)  measurements.  The  altitude  dependence  of  the  N(^D)  column  density 
measurement,  which  looks  18®  into  the  afterglow,  is  controlled  by  quenching  of  N(^D)  by 
secondary  electrons  with  increasing  importance  of  the  N(^D)+02  reaction  near  100  km.  On 
upleg,  the  largest  discrepancy  between  the  model  and  data  is  at  109  and  111  km  on  upleg,  which 
is  not  understood,  although  it  is  likely  related  to  the  uncertainty  in  the  secondary  electron 
quenching  rate  of  N(^D).  The  differences  between  upleg  and  downleg  N(^D)  column  densities 
simply  reflect  the  differences  the  ion  pair  production  rate,  although  the  ratio  of  upleg  to  downleg 
N(^D)  column  densities  is  less  sensitive  to  ion  pair  production  rate  due  to  the  significant 
quenching  of  N(^D)  and  the  fact  that  the  primary  production  process  for  N(^D)  is  recombination 
of  secondary  electrons  with  Nj.  The  production  rate  of  N(^D)  from  recombination  roughly 
increases  as  the  square  root  of  the  ion  pair  production  rate. 

The  comparison  of  the  model  predictions  with  the  3466  A  data  shown  in  Figure  8  indicate 
that  the  N(^P)  chemistry  is  not  as  well  established  as  the  N(^D)  chemistry.  This  comparison 
would  indicate  that  the  assumed  production  rate  of  0.22  N(^P)/ion  pair  used  in  Reaction  (4b)  is 
approximately  a  factor  of  2  low.  Since  the  total  production  rate  of  N  atoms  is  fairly  well 
established,  increasing  the  rate  for  Reaction  (4b)  by  a  factor  of  2  would  imply  that  only  N(^P) 
is  produced  by  dissociation  of  N2  by  primary  electrons.  Such  a  modification  in  the 
N(^D)/N(^P)  relative  production  rate  from  Reaction  (4)  would  not  significantly  alter  the 
excellent  agreement  between  the  predictions  and  the  5200  A  data  since  the  major  production 
mechanism  for  N(^D)  is  the  recombination  of  Nj  with  secondary  electrons.  Reaction  (6),  as 
shown  in  Figure  1.  Furthermore,  the  increase  in  N(^P)  would  provide  more  N(^D)  due  to  the 
increased  quenching  rate  from  Reaction  (14). 
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Figure  5. 


Figure  6. 


Total  Number  of  Ion  Pairs  Along  the  Line-of-Sight  as  a  Function  of  Mission 
Elapsed  Time  for  the  EXCEDE  III  Experiment. 


Calculated  Dose  Time  as  a  Function  of  Mission  Elapsed  Time  for  the  EXCEDE  III 
Experiment. 
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Figure  7,  Comparison  of  the  Time  Dependence  of  the  Calculated  N(^D)  Column  Density  (o)  with 
the  5200  A  Photometer  Data  (•). 
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Figure  8.  Comparison  of  the  Calculated  Time  Dependence  of  the  N(^P)  Column  Density  (o)  with 
the  3466  A  Spectrometer  Data  (•). 
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5.  SUMMARY 


The  EXCEDE  III  experiment  has  provided  quantitative  characterization  of  the  N(^D)  5200 
A  and  N(^P)  3466  A  emissions  resulting  from  the  irradiation  of  the  lower  thermosphere  (90-115 
km)  by  an  intense  electron  beam.  A  detailed  chemical  kinetics  model  indicates  that  the  N(^D) 
emission  can  be  quantitatively  explained  by  assuming  accepted  production  rates  of  N(^D)  and 
N2  created  by  the  electron  beam  dissociation,  consistent  with  previous  auroral  models.  Using 
the  previous  estimate  of  0.22  N(^P)  per  ion  pair^^'*^  from  electron  dissociation  of  N2  results  in 
a  model  prediction  which  is  approximately  50%  below  the  data,  in  disagreement  with  previous 
analyses  of  N(^P)  auroral  and  dayglow  emission.^^^’^^’^^^  However,  it  should  be  recalled  that 
the  EXCEDE  III  measurement  of  N(^P)  was  not  made  under  steady  state  conditions,  which  are 
appropriate  to  the  atmospheric  observations.  The  present  analysis  would  imply  that  the  amount 
of  N(^P)  from  N2  dissociation  should  be  about  a  factor  of  2  greater  than  previously  assumed  or 
an  additional  source  of  N(^P)  exists  for  intense  electron  beams  such  as  that  used  in  EXCEDE 
III. 
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